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Abstract. We report BeppoSAX observations of spatially re- 
solved hard X-ray emission from IC 443, a supernova rem- 
nant interacting with a molecular cloud. The emission is shown 
to come from two localized features spatially correlated with 
bright molecular emission regions. Both hard X-ray features 
have soft X-ray counterparts, in one case shifted by ~ 2' toward 
the remnant interior. The spectra of X-ray emission from both 
isolated features have photon index < 2 in the MECS regime. 
The emission detected from the remnant with PDS detector ex- 
tends up to 100 keV. We discuss the observed properties of the 
hard X-ray features in relation to non-thermal emission from 
shocked molecular clumps and pulsar wind nebula. 

Key words: Acceleration of particles; Radiation mechanisms: 
non-thermal; ISM: clouds; ISM: cosmic rays; ISM: individual 
object: IC443; ISM: supernova remnants 



1. Introduction 

Massive stars that are the likely progenitors of core collapsed 
supernovae are expected to be spatially correlated with molec- 
ular clouds. IC 443 (G189. 1+3.0) is one of the best space labo- 
ratories to study rich phenomena accompanying supernova in- 
teraction with a molecular cloud. This was established from 
observed shock-excited molecular line emission from OH, CO 
and H 2 (e.g. Burton et al. 1988; van Dishoeck et al. 1993; Ce- 
sarsky et al. 1999). Observations of OH masers and CO clumps 
towards IC 443, W28 and W44 also provide convincing evi- 
dence of molecular cloud interaction in this case (Claussen et 
al 1997). 

IC 443 was a target of X-ray observations with HEAO 1 
(Petre et al. 1988), Ginga (Wang etal. 1992),,/?<9&4r(Asaoka 
& Aschenbach 1994) and ASCA (Keohane et al. 1997, K97 
hereafter). The soft X-ray 0.2-3.1 keV surface brightness map 
of IC 443 from the Einstein Observatory (Petre et al. 1988) 
shows bright features in the northeastern part of the remnant as 
well as bright soft emission from the source central part. The 
presence of nearly uniform X-ray emission from the central 
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part of the remnant was clearly seen also by ROSAT (Asaoka & 
Aschenbach, 1994), and corresponds to the emission from hot 
(T~ 10 7 K), low density gas interior to the shock. ASCA GIS 
observations by K97 discovered the localized character of the 
hard X-ray emission. They concluded that most of the 2-10 keV 
GIS photons came from an isolated emitting feature and from 
the southeast elongated ridge of hard emission. Preite-Martinez 
et al. (1999) have reported a hard (> 14 keV) component with 
BeppoSAX/PDS and two hot spots corresponding to the ASCA 
sources of K97 using the MECS. 

The supernova remnants (SNRs) that are likely candidates 
to be 7-ray sources in CGRO EGRET observations (Esposito 
et al. 1996) also show evidence for interaction with molecular 
gas. In the case of IC 443, the molecular line emission region 
is partially inside the EGRET 7-ray detection circle. 

Theoretical models of SNRs interacting with clouds were 
discussed recently by Sturner et al. (1997), Baring et al. 
(1999), Chevalier (1999). The non-thermal multi-wavelength 
spectrum of a SNR interacting with a molecular cloud was 
studied by Bykov et al. (2000). They showed that the propa- 
gation of a radiative shock wave within a molecular cloud lead 
to a substantial non-thermal emission both in hard X-rays and 
in 7-ray. The complex structure of molecular cloud consisting 
of dense massive clumps embedded in the inter-clump medium 
could result in localized sources of hard X-ray emission corre- 
lated with both bright molecular emission and extended source 
of non-thermal radio and 7-ray emission. In this Letter, we 
present the archive BeppoSAX observations of IC443 and, in 
particularly, we discuss the features detected in the MECS hard 
X-ray mosaiced map of the remnant. 

2. Observations and data analysis 

Results from the the Medium-Energy Concentrator Spectrom- 
eter (MECS; 1.8-10 keV; Boella et al. 1997) and the Phoswich 
Detection System (PDS; 15-300 keV; Frontera et al. 1997) on- 
board BeppoSAX are presented. 

We have used all the 4 observation of IC443 available in 
the public BeppoSAX archive; Table |l] lists their coordinates 
and exposure times. Good data were selected from intervals 
when the elevation angle above the Earth's limb was >4° and 
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Fig. 1. MECS vignetted and exposure corrected mosaic of the four IC443 observations listed in Table |T| Grayscale is linear and 
contours are 20%, 40%, 60% and 80% of the peak (3.6 and 2.2 x 10 -3 cnt s _1 arcmin -1 for the soft and hard image, respectively). 
In the right panel, the sources discussed in the text are marked A, B and C, and dashed contour correspond to a v — 1 — 05(1) 
molecular H line intensity of 6.6 x 10~ 14 erg cm~ 2 s _1 reported by Burton et al. (1988); arrows indicate the location of peaks 
with an H line intensity at least 14 times the contour levels. 



Table 1. Public BeppoSAX Observations of IC443 



Obs. 


1 exp 

ksec 




Coord (J2000) 




Date 


IC443H 


21.5 


6' 1 


17 m 05 s ±22 d 22 m 


19 s 


11/4/1999 


IC443E 


48.4 


6" 


18 m 00 s ±22 d 34 m 


01 s 


1 1/4/1998 


IC443C 


67.1 


6" 


16 m 01 s ±22 d 28 m 


55 s 


18/10/1997 


IC443N 


23.8 


6 h 


17 m 12 s ±22 d 48 m 


50 a 


20/10/1997 



when the instrument configurations were nominal, using the 
SAXDAS 2.0.0 data analysis package. 

Fig. [] shows the MECS mosaics of IC443 in different 
bands. The images are vignetting and exposure corrected. The 
2.0-4.0 MECS images show an emission morphology similar 
to the GIS broad-band maps presented by K97. On the other 
hand, the hard MECS image above 4 keV is very different. 
There is weak diffuse extended emission in the north of the 
remnant, and two strong compact sources, Src A and B. Table 
U reviews the properties of the two sources. In Fig. [lj we also 
marked the weak extended hard X-ray emission with C. Both 
the Src A and B were observed by K97; Src A corresponds to 
the hard source named "HXF", while Src B corresponds to the 
hard extended emission called "The Ridge" (also reobserved 
by Olbert et al. 2000, and detected by Preite-Martinez et al. 
1999 as 1SAX J0618.1+2227). In the hard X-ray image in Fig. 
[j] we have also reported the contour of the molecular hydrogen 
line emissivity observed by Burton et al. (1988), which shows 
that the hard X-ray sources are closely correlated with it. 

In Fig. ||we report the 4.0-10.5 keV profile of Src A and 
Src B compared with the MECS PSF profile in the same energy 



Table 2. Compact sources detected in the hard MECS images. 



Name 


Count-rate[| 


Position 








(s" 1 ) 


(J2000) 






Src A 


3.7 ±0.2 x 10~ 2 


6" 17 m 6.3 s ±22 d 


21 m 


37 s 


SrcB 


5.9 ±0.6 x 10" 3 


gh lg m g 7 s +22 d 


27™ 


28 s 



In the 4.0-10.0 keV energy range. 



Table 3. Summary of MECS spectral fitting results. 



Model 


kT 

keV 


7 


X z /dof 






Src A 




PL. only 




1 96+ 021 


34/41 






SrcB 




PL.+MEKAL 


1MI -0.31 


1 c+0.9 
l -°-0.9 


41/33 



range. The PSF profile was derived from a BeppoSAX obser- 
vation of the very bright point source Cyg X- 1 . Src A seems to 
be extended, while Src B is consistent with the expected PSF. 
The latter was reported as an elongated ridge by K97, but we 
note that the ASCA GIS pointing direction is such that the loca- 
tion of Src B is seen at large off-axis angle, which can account 
for most of source extension observed by ASCA. Instead, the 
BeppoSAX Src A and B are placed only at 5' and 1' off-axis, 
where the instrumental distorsion is at its minimum. 

We have also investigated if Src A and Src B have soft X- 
ray counterpart using archive ROSAT observations. Fig. |]is a 
PSPC 0.2-2.0 keV image of the region containing MECS Src 
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Fig. 2. Hard X-ray profile (4.0-10.5 keV) of the MECS Src A and B (dashed line) compared with the expected profile of point-like 
source (solid line) derived from a long MECS observation of the X-ray binary Cyg X-l. 



A and B, with the MECS 4.0-10.5 keV contours of Fig. [I] (right 
panel) overimposed. Src A has an HRI counterpart quoted by 
K97 and the PSPC WGACAT source 1WGA J0617. 1+2221 
with a rate of2.8±0.3x 10~ 2 cnt s _1 is located inside the 
BeppoSAX error circle. Src B has a weaker PSPC counterpart 
indicated by the arrow in Fig. |], located ~ 1.5' north of the 
MECS source centroid and not reported in other PSPC cata- 
logs. We have verified that the extended hard X-ray emission 
we have labelled with C in Fig |l] has a thermal origin and it 
represents the weak hard tail of the bright soft X-ray nebula 
visible in the PSPC image. 

We have performed spectral fits of the MECS spectra of Src 
A and B using the MEKAL optically thin thermal plasma model 
modified by the interstellar absorption assuming the standard 
abundances. Spectra were extracted from circular regions with 
radii of 4' and 3' for Src A and B, respectively. We have used 
the IC443H observation for Src A and IC443E observation for 
Src B. Spectra have been rebinned to have at least 20 counts 
per channel. The MECS background was collected in an annu- 
lus between 5' and 8' for Src A, while we have used the stan- 
dard background for Src B because no source-free background 
region could be found around it. 

A thermal origin for the spectrum of Src A is strongly re- 
jected (x 2 /d.o.f. = 111/41), while a power law model nicely 
fit the data (Table |3|). The spectrum of Src B cannot be repro- 
duced by a power law or by a thermal emission model alone 
(X 2 /d.o.f. =61/33 and 78/33 respectively) but requires a com- 
bination of the two, reported in Table [3| 

We have also analyzed the PDS spectra. Background was 
subtracted using the OFF collimator positions; we have ver- 
ified that no other contaminating sources are present in the 
background positions. Spectra were collected using the Vari- 
able Rise Time Threshold mode and have been rebinned fol- 
lowing a logarithmic scheme suggested by the PDS hardware 
group to not undersample the instrument spectral resolution. 
Table Q reports the observed countrates and flux in different en- 
ergy ranges. A joint fit of Src A MECS spectrum and IC443H 
PDS spectrum (shown in Fig. ||) yields best-fit results similar 



MECS and PDS Spectrum of Source A 
Power— law mode 




10 100 
Energy (keV) 



Fig. 3. MECS and PDS spectrum of Src A. The power law 
model which best fit the data (Table ||) is also reported with 
residuals. 

to the ones reported in Table [| and a PDS/MECS normalization 
ratio of 1.40 ± 0.45. This value is consistent with the expected 
value of 1 , also considering that the expected constribution of 
Src B in the PDS spectrum of IC443 could be up to ~ 50%. As 
for Src B, the joint fit yields a PDS/MECS ratio of 2.8 ± 2.5, 
we argued that the PDS spectra of IC443E is heavily contam- 
inated by the brighter Src A. In fact, using the best-fit model 
of Src A and considering the triangular spatial response of the 
PDS collimator, we predict a Src A contribution of ~ 6 x 10 -2 
cnt s _1 in the PDS spectra of Src B between 15 and 30 keV, 
consistent with the observed rate (Table Therefore, a proper 
analysis of the PDS spectra of Src B cannot be done. 

3. Discussion 

The compact hard X-ray features may be interpreted in terms 
of the interaction between the remnant shock and a molec- 
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Table 4. Observed background subtracted count-rates and un- 
absorbed fluxes in different energy ranges, MECS for E < 10 
keV and PDS for E > 10 keV. 



Band 


Count-rate 


Flux xl0~ 12 (fl) 


keV 


(10 -2 s" 1 ) 


erg cm -2 s~ 






Src A 


2-10 


8.2 ±0.2 


7.5 ±1.2 


15-30 


11.5 ±3.0 


9.1 ±2.4 


30-100 


7.2 ±4.3 


11.1 ±6.6 


100-220 


5.8 ±3.7 


41.5 ±26.7 






SrcB 


2-10 


2.0 ±0.1 


1.8 ±0.20 


15-30 


4.3 ±2.4 


3.4 ± 1.9 


30-100 


< 3.9 


< 1.1 



8 Computed using the best-fit models and 
values reported in Table [| 
* 47% MEKAL, 53% P.L. 



ular cloud. Bykov et al. (2000) have presented a model of 
non-thermal emission from an evolved SNR interacting with 
a molecular cloud, but the predicted diffuse hard X-ray emis- 
sion in the MECS bandwidth is below the sensitivity threshold. 
On the other hand, several clumps have been observed in the 
molecular cloud near IC443. The spatial scales of the molec- 
ular clumps from IR line observations are in the range from 
150" down to 1". The non-thermal X-rays from clumps were 
modelled by Bykov et al. (2000), who predicted a hard spec- 
trum with a photon index below 2.0. The emission is originated 
from electrons of energy below 1 GeV. The predicted flux den- 
sity is - 6 x 10~ 5 keV cm~ 2 s" 1 keV" 1 at 10 keV for a 30 
km s~ 4 shock in a clump of half parsec radius and density of 
10 4 cm -3 . Using a power law with 7 = 1.5, this corresponds 
to - 3 x 10~ 2 MECS cnt s" 1 in the 4.0-10.5 keV energy 
range and ~ 10~ 2 PDS cnt s" 1 in the 15-30 keV. The pre- 
dicted count-rates are consistent with the observed rate of Src 
B (Table |J), having in mind that 10 4 cm~ 3 density is a rep- 
resentative number and other parameters could easily account 
for the difference. Because of hard spectrum and heavy absorp- 
tion, localized spots of few arcmin size would be seen only in 
hard X-rays. Thus, the compact Src B correlated with bright 
spot of molecular hydrogen emission (Fig. [j]) can be shocked 
molecular clump. The observed soft X-ray brightening having 
apparent shift from hard X-ray features towards the SNR in- 
terior (Fig. Q), can be attributed to the effects of the shocked 
clump edges. The morphology and the spectrum of the Src A 
are somewhat different from that of Src B. 

As already discussed by K97, a low-luminosity pulsar neb- 
ula is a plausible explanation for isolated hard X-ray sources 
(especially for the slightly extended Src A), but we do not ex- 
pect to see two of them in one snrQ. High resolution Chandra 
or XMM data are needed to distinguish between pulsar neb- 

1 However, Asaoka & Aschenbach (1994) pointed out that the soft 
X-ray emission is consistent with the presence of two large SNRs, 10 
and 10 5 yr old. 
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Fig. 4. 0.2-2.0 keV PSPC image of the region cotaining the two 
MECS hard X-ray sources. The hard X-ray contours of Fig. [j] 
are overlaid. The arrow shows the soft X-ray counterpart of Src 
B. 

ula and shocked clump interpretation. We have also performed 
a timing analysis of MECS events collected inside a circle of 
3' radius in the 4.0-10.5 energy range for both Src A and B, 
and we found no pulsations at 99% confidence level, with an 
upper limit of the pulsed fraction of a sinusoidal signal in the 
10~ 3 — 64 Hz frequency range of 5% and 6%, respectively, 
more stringent than the one derived by K97 using ROS AT data. 

We have also estimated the Src A and B EGRET 7-ray E > 
100 MeV flux contribution, using the best-fit models reported 
in Table | and they are 3.0± ;| x 10~ 8 and lOOtjjf x 10~ 8 
photons cm~ 2 s . The expected 7-ray flux of Src B therefore 
is compatible with the value of 50.0±3.0 x 10 -8 photons cm -2 
s _1 observed by Esposito et al. (1996). Gamma-ray observa- 
tions with forthcoming INTEGRAL mission with expecting an- 
gular resolution about 12' could better constrain the spectrum 
of the sources above 30 keV. 
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